Solitons of the Resonant Nonlinear Schrodinger 
Equation with Nontrivial Boundary Conditions 
and Hirota Bilinear Method 



Jyh-Hao Lee 1 and Oktay K. Pashaev 2 
1 Institute of Mathematics, Academia Sinica 
Taipei, 11529, Taiwan 
2 Department of Mathematics, Izmir Institute of Technology 
Urla-Izmir, 35430, Turkey 

February 5, 2008 



Abstract 

Physically relevant soliton solutions of the resonant nonlinear Schrodinger 
(RNLS) equation with nontrivial boundary conditions, recently proposed 
for description of uniaxial waves in a cold collisionless plasma, are consid- 
ered in the Hirota bilinear approach. By the Madelung representation, the 
model transformed to the reaction-diffusion analog of the NLS equation 
for which the bilinear representation, soliton solutions and their mutual 
interactions are studied. 



1 Introduction 

Recently for description of low dimensional gravity (the Jackiw-Teitelboim model), 
and a medium response to the action of a quasimonochromatic wave with com- 
plex amplitude ip(x,t), which is slowly varying function of the coordinate and 
the time, a novel integrable version of NLS equation namely 

was introduced pQ. This has been termed the resonant nonlinear Schrodinger 
(RNLS) equation. It can be considered as a third version of the NLS equation, 
intermediate between the defocusing and focusing cases. Even though the RNLS 
model is integrable for arbitrary values of the coefficient s, the critical value 
s = 1 separates two distinct regions of behavior. Thus, for s < 1 the model 
is reducible to the conventional NLS equation. However, for s > 1 it is not 
reducible to the usual NLS equation, but rather to a reaction-diffusion (RD) 
system. In this case the model exhibits resonance solitonic phenomena pQ. 



The RNLS equation can be interpreted as a particular realization of the NLS 
soliton propagating in the so called "quantum potential" Uq(x) — IV'Ux/IV'l- 
This potential, responsible for producing the quantum behavior, was introduced 
by de Broglie [2] and was subsequently used by Bohm [3j to develop a hidden- 
variable theory in quantum mechanics. It also appears in stochastic mechanics 
[4 . Connections between such non-classical motions with the internal spin mo- 
tion and the zitterbewegung have been considered in a series of papers (see [5]). 
Quantum potentials also appear in proposed nonlinear extensions of quantum 
mechanics with regard both to stochastic quantization [B] , [7] and to corrections 
from quantum gravity [5]. It is noted that the RNLS equation, like the conven- 
tional NLS equation, may also be derived in the context of capillarity models 

m, ma 

Very recently [11] it was shown that the RNLS equation appears in plasma 
physics, where it describes the propagation of one-dimensional long magnetoa- 
coustic waves in a cold collisionless plasma subject to a transverse magnetic 
field. The complex wave function satisfying RNLS equation is a combination of 
plasma density and the velocity fields as in the Madelung representation. The 
Backlund-Darboux transformations along with novel associated nonlinear su- 
perposition principle were presented and used to generate solutions descriptive 
of the interaction of solitonic magnetoacoustic waves. This application requires 
to consider a solution of the RNLS equation with nontrivial boundary condi- 
tions at infinity. The goal of the present paper is to derive such solutions in the 
Hirota bilinear approach and study their mutual interactions. 



2 Magnetoacoustic waves in cold plasma 

The dynamics of two-component cold collisionless plasma in the presence of an 
external magnetic field B [12] , [13] for uni-axial plasma propagation 

u = u(x,t)e x , B = B(x,t)e z (2) 

reduces to the form [14] 

5 + sW"* ,3) 



du du B dB 
u 



= 0, (4) 
B-p. (5) 



dt dx p dx 

d_ fldB y 
dx \p dx 

(where we set B — 1 and plasma density p = 1 at infinity) This system is 
equivalent to that of Whitham and has also been derived by Gurevich and 
Meshcherkin [T5]. It describes the propagation of nonlinear magnetoacoustic 
waves in a cold plasma with a transverse magnetic field. It has been shown 
recently by El, Khodorovskii and Tyurina [19] that a system of the type ©-(O 
also occurs in the context of hypersonic flow past slender bodies. 
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3 A shallow water approximation 



Here, we consider a shallow water approximation to the magnetoacoustic system 
([3])-([5]). Thus, rescaling the space and time variables via x' = fix and t' = fit, 
we have 

BOB 
p dx' 

d (\dB" 



du du 
dF +U dx 7 







B-p 



dx' \ p dx' t 

On expansion of B as a power series in the parameter f3 2 according to 

B = p + [l 2 b 2 ( Pl p x ,,p x , x ,,...) + 0(p 4 ), 
insertion into ijHJ yields 



dx' \ p dx' 



Substitution of © into (J7J) yields 



du du dp 2 

dF + u M + w +(3 



1 d 3 p 2 dp d 2 p 
p dx' 3 p 2 dx' dx' 2 



1 dp_ 

p dx' 



= 



to O{0 2 ). Accordingly, the following system results: 

dp ® f \ 



du du 
dt i + U dx 7 



dt' dx' 

^ + (3 2 — 
dx' dx' 



1 d 2 p I (1 dp"" 2 



p dx' 2 2 \p dx' 



= 0. 



(6) 
(7) 
(8) 



(9) 
(10) 

(11) 

(12) 
(13) 



This describes the propagation of long magnetoacoustic waves in a cold plasma 
of density p moving with velocity u across the magnetic field as given by @, 
©(see [TB], [IT]). In this system the dispersion is negative, i.e., the wave velocity 
decreases with increasing wave vector k. 



4 Resonant NLS 

Introducing the velocity potential S(x,t) = — \ J x u(x, t)dx so that u 
and integrating once equation (fT3"|) we get the system 



8S 



dS 



-m + {dx-) + 2 P + Y 



dp 

dt 
1 



2 — (p—) 

dx dx 







1 d 2 P 1/1 dp 
p dx 2 2 \p dx 



-2dS/dx 

(14) 
(15) 
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(here we skip prime superscript). Combining p and S as one complex function 

V = Vpe-' iS (16) 

the system (IT41 . (fT5j) can be represented as the nonlinear Schrodinger equation 
with quantum potential, the RNLS dXJ) , 

where A = —2, s = 1 + /3 2 , and we have used expression for the quantum 
potential 

1 d 2 p 1 (I d P \ 2 1 d 2 Jp 



p dx 2 2 \p dx) yfp dx 2 ^ ^ 

Since parameter s > 1 equation (fT7|) can not be transformed to the NLS pQ. 
But by combining p and 5 as a couple of real functions 

e(+) = \/pe^ S , e M = ~VP e ~^ S ( 19 ) 

so that e^ + ' > 0, e^ - ) < 0, equations (fl"4f . (fl"5)) can be written as the reaction- 
diffusion (RD) system 

where r = fit' . 

Linearization of JT7]) near the "condensate" solution ip — gives 
dispersion lo = ^fpok — |^ k 2 . This dispersion is negative, i.e. the wave 
velocity decreases with increasing wave vector k and is unstable for short waves 
with k > k cr , where k cr — ./po//3. However this instability results from the 
truncation of the dispersion relation 

Po 

for the system ©, Q, © and does not correspond to any actual physical effect 
for shallow water waves [18J . Though this system is linearly stable for all wave 
numbers k, it is not known to be integrable. Therefore, it is not as suitable 
for studying wave interactions as the RNLS, which is completely integrable and 
admits a rich variety of the exact solutions. 



5 Steady State Flow and Solitons 

Now we consider the system (fT2"|) . (|13[) for the steady state flow. It describes 
the motion with fixed velocity u(t, x) = uq = const for which the continuity 
equation (|12|) implies + ^off — or that fluid density has the travelling 
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wave form p — p(x — uot), where x' — f3x, tf — (3t. Equation (|13[) in this case 
gives 



1 d 2 P 1/1 dp s 1 



p dx 2 2 \p dx 



= a = const. 



(22) 



This has a simple physical interpretation. The motion with fixed velocity implies 
that sum of all forces acting on the system is zero. In our case Eq. (j22"]) describes 
compensation of the nonlincarity by the quantum potential. Using (|18p this 
equation gives 



yfp dx 2 

In terms of y(z) = y/~p, z = x — u t, then we have nonlinear equation 

d 2 v a 1 , 

or multiplying with y' and integrating once 

2 



dz 

This equation has solution 



1 



-y = 2b = const. 



y = 2pdn[p(x — uot), k] 



(23) 



(24) 



(25) 



(26) 



where dn is the Jacobian elliptic function with the modulus k, and p is an 
arbitrary constant. It is connected with integration constants a > and b < 
by equations 

w = l+ *°- w—" (27) 

where k' = \J\ — k 2 is the complimentary modulus of the Jacobian elliptic 
function. These equations fix relation between the constants 



a = 2p z 1 - 



2p 4 



(28) 



and the modulus k = yl + b/2p 4 . Then for density p we have the travelling 
wave solution 

p{x,t) — Ap 2 dn 2 [p(x - uot), k] . (29) 

With fixed potential from equation JTSJ) and (??) for u Q = -(2/ 0)dS/dx, 
we have the Hamilton- Jacobi equation 



/3 dt ' (3 2 
which has solution in the form 

»o 



1 dS 1 fdS . 



dx 



S = [3 



So 



(30) 



(31) 
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or using 



S = f3 



S -^x + (^ +P 2 (2-n 2 ))t 



(32) 



For the degenerate case k — 1 which corresponds to b — and as follow a = 2p 2 , 
the elliptic solution (|2^|) . (f3"2")) becomes the soliton 



p(x,t) = 4p 2 secft. 2 [p(a; — Uoi)] 



S = /9 



c u 



P 2 )t 



(33) 
(34) 



More general form of the travelling wave appears if we consider solution of 
the system (1121) . (113|) in the form u(x, t) =u(x — Uot), p(x,t) = p(x — Uot). Then 
the first equation (fl"2"|) implies 4z\{u — Uo)p] = 0, where z = x — u t and 



u = Uq + 



C 



(35) 



where C = const. Substituting to the second equation (flU)) and integrating once 
we have 



1C 2 
2 p A 



i d 2 p i n dp 

p dz 2 2 \p dz 



A = const. 



(36) 



Multiplying with p 2 and differentiating once we have (travelling wave form of 
the KdV equation) 



'd 3 p dp 0A dp\ 



which implies after one integration 
d 2 p , 3 2 



This gives equation 



and solution 



dz 2 



2Ap = B = const. 



-p 3 + 2Ap 2 + 2Bp + C 2 



(37) 



(38) 



(39) 



p(x — Ugt) = ax + (0*3 — ai) dr7, 2 [-(a 3 — ai) 1 ^ 2 (x — u t), n] (40) 

with modulus of elliptic function k 2 — («3 — a^) /(a3 — ai), constants ai, ct2, a 3 
and 

(ai a 2 a 3 ) 1/2 



u(x — Uot) = Uq 



P 



(41) 



This solution was first reported in [17] • We have particular reductions: 



G 



1) at\ = aa = and corresponding n 2 = 1 so that (|4*tJ)) reduces to (f^^j) . 
with following identification of parameters 0*3 — 4k 2 , b — —a2k 2 /2 and velocity 

U = Uq. 

2) ai = a2 7^ 0, then again k 2 = 1 but solution 

p(x - u t) = ai + (a 3 - ax) sech 2 [^(a 3 - ai) 1/2 (x - u t)\ (42) 

has nontrivial asymptotics. The physical value relevant for plasma physics is 
a% = 1 leading to limui^oo p = 1. Denoting a 3 = a 2 we can write it in the form 

p(x - uot) = 1 + ^= - (43) 

PK ' cosh 2 [^3^ - uot)] V ' 

Using the above results we can now construct solutions of RNLS (fl"7|) . Sub- 
stituting Eqs. (f29|)(f32"l) to Eq. (fT6)) and changing parameter k = p//3, we have 
quasi-periodic solution 

t') = 2[3k dn[k(x' - u t'), K ] e ~^o-^x'+i4+P 2 k 2 (2- K 2 )]t'} (44) 

where x' = f3x, t' = (it. In the limit k = 1 it gives the envelope soliton solution 

- l {4, -^. x ' + [^+f3 2 k 2 ]t'} 

W) = W eosh^-^0 (45) 

For the reaction-diffusion system ([20]) correspondingly, we have the " dissipative" 
periodic solution 

e i±} {x',T) = ±2/3k dn[k(x' - vt), K ] e M4>o-^x'+[4+kH2-^)]r} ^ 

where velocity v = uq//3, k = p/j3 and the dissipative analog of the envelope 
soliton 

e ±{<Po-^x' + [^-+k 2 ]r} 

e^(x',r) = ±2(3k- — — — (47) 

coshfc(x' — VT) 

the so called "dissipaton" solution [1] 



6 Bilinear form and solitons. Trivial Boundary 
Conditions 

The reaction-diffusion representation (|2U)) of the RNLS (fTT|) is the key point for 
construction multisoliton solutions. Due to algebraic similarity of RD system 
with NLS it is easy to write the bilinear representation for it. Representing two 
real functions e^ + \ e^ - ' in terms of three real functions 



c 



(±) 



2/3^ (48) 
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we have the next bilinear system of equations 

(±D T - D 2 X ){G {±) ■ F) = 0, D 2 X (F ■ F) = -2G {+) G { - ] (49) 
The corresponding solution of the RNLS (jTTJ) has 

Mx, t)| 2 = P = -eWe(-) = 2f^H = 4/3 2 ^ (50) 



The one-dissipaton is given by the following solution of system ([49]) : G 



±e v i, F = 1 + eX +J 'i e^ 1 ' 1 = (fc+ + /cf)- 2 , where ^ = jfefz ± (fc±) 2 T + 

TJi and fcf, t?^ are constants. In terms of redefined parameters, k = (k^ + 
ki)/2, v = —(kf — fcj") it acquires the form l[47p. In the space of parameters 
(u, fc) there exist the critical value v cr u = 2k for solution (??) so that when 
v < v cr it, one has e — > at infinities, so the vanishing b.c. for the dissipaton. 
At the critical value the solution is a kink steady state in the moving frame 

= ±ke ±k £°(l =F tanh fc£), with constant asymptotics e ± — ► ± 2ke ±k ^° for 
a; — * =Foo and e — ► ± for x — > ±00. In the over-critical case v > v cr u, 

— * ±00 for x — ► =Foo and e 1 * 1 — * ± for x — > ±00. 

For the two-dissipaton solution we have 

G ± = ±[ e vt + e <£ + ( ^ + _ f^t+^+nf + ( ± fc | 2± + _ ) 2 e ^+^+^] ; (51) 



F = 1 



e^+ni e nt+v 2 e vt+V! e v}+ii 2 



(fc n ) 2 (fc 12 ) 2 (fc 21 ) 2 (fc 22 



-^2 



"12 h, 2i 



+ 



where fcf/ = fc? + fc), fcf/ = fcf-fcj, Tjf = kfx±(kf) 2 T + r] ± ^ . This solution 
shows the resonance character of dissipatons interaction pQ 

7 Bilinear form and solitons. Non-trivial Bound- 
ary Conditions 

As was first noticed by Hirota, for NLS equation of defocusing type with nonva- 
nishing boundary conditions, the bilinear form of equations should be modified 
[20] , After substituting representation ([48]) to the system ([20]) the decoupling 
system is choosen in the form 

(±D T -D 2 X + A)(G ± • F) = 0, {Dl - A)(F • F) = -2G+G~ (53) 



where we have introduced constant A to be determined. Equation (|48|) and the 
second one of ([53]) imply — e^eS^ — —4f3 2 [-| — (lnF) xa; ]. Expanding G ± and 
F in Hirota's power series 

G ± = ±gf{l + eg? + ^dt + .-.), F = l + e/i + e 2 / 2 + (54) 
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and requiring lim| a; |^ 00 (ln.F) xa: = Owe have the boundary condition 

A 



ai = lira [-e (+) e (_) ] = lira -4/3 2 

\x\— >oo |cc|— >oo 



2 (inn 



2(3 2 X (55) 



which fixes constant A = — oti/2/3 2 . In the zero order approximation we have 
the system 

(±D T -Dl + A)(.g ± • 1) = 0, (D 2 X - A)(l • 1) = 2g+ g^ (56) 

It has a solution in the form g^ = /3 ± e e± , where 0q = (3 + (3~ = —A/2 = ai/4/3 2 , 
± = ±fcx ± (fc 2 — A)t, /J* = /3oe ±7 ° . Using following properties of Hirota's 
derivatives 

D x (.fg-h) = ^-gh + fD x (g-h) (57) 
f) 2 f f)f 

D l(f9 ■ h) = g^gh + 2-±D x {g ■ h) + fD 2 x (g ■ h) (58) 
the bilinear system becomes 
( T D T ± 2kD x + D 2 X ){{1 + eg? + e 2 gf + ■■■) ■ (1 + e/i + * 2 / 2 + •••)) = 0, (59) 
( J D 2 +2/3 2 )((l + e/ 1 + ...)-(l + e/i + ...))-2/3 2 (l + e.g+ + ...)(l + e.gr + -)- (60) 

7.1 One soliton solution 

In the first order we have the system 

( T d T ± 2fc9 x + d 2 x )g± + (±0 T T 2kd x + d 2 ).fr = 0, (61) 

{d 2 x + 2(3t)h=Pl{gt+g- l ). (62) 

Considering solution in the form gf = a^e m , fi — bie' ni , where 771 = k\x + 
uj y t + rfl we get = 7161, a~ = j^b 1 , 71 = (wi - 2kk 1 + fc 2 )/(wi - 2kk x - fc 2 ), 
where dispersion formula is 

ut = ki (2k ± ^k 2 +Al3 2 ) (63) 

It is worth to notice that in contrast to the dark soliton of the defocusing NLS 
equation, in our case no restrictions on values of k\ appear. 

Truncation of the Hirota expansion at this level gives one dissipative soliton 

e (±) = ±203 ± c ±lkx+(k2+2f32n)T] 1 + 7l±le??1 (64) 

l + e m K J 

were we have absorbed constant b\ to the exponential form f\\ = k\x + uj\t + 
rfl + In bi . Then we have one soliton density 
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This solution can be represented in the form 



e (±) =± x /SI/i ±1 e 



±1 ±[kx+(k 2 +2/3^)r] 



7 ±x + l | 7 



±1 



1 77 

tanh — 

2 2 



or 



e (+) - +V^T £ e +[^+(fe 2 + 2 ^ 2 )-] ( 7 + 1 + ( 7 - 1) tanh 



e (") = _ /H7A e -[^+(fc 2 +2/3 2 W 



and the product is 



2p, 



p = — e^ + M ' = a\ 



11 fi 
- + 1 + ( l)tanh- 

7 7 2 



1 + 



(7 -If 
47 cosh 2 § 



with asymptotic limi a .i_ >0O p —* ot\. Explicitly it is 



p = — e^ + M ' = ai 



1 + ^ ^ 2 | (z + (2k ± V /fc 1 +4/? 2 )r + 



.TO 



where /3 2 = l/(4/3 2 ai). For the velocity field we have 



P (-) P (+) 



U = 



= -2/c- 



(7 2 - l)*i 



e(-) e(+) " (7 -l) 2 +4 7 cosh 2 f 

Let us consider particular solution for k — 0, then as follows 



wi = ±fci^/fc 2 +4/3, 



(66) 

(67) 
(68) 

(69) 

(70) 

(71) 
(72) 



is the Bogolubov dispersion from the theory of superfluidity of a weakly non- 
ideal Bose gas. For k\ » 1 it is of the non-relativistic free particles form 
0J1 i=a kf, while for k\ « 1 it is of the relativistic collective form lo\ rs 2/3 ki. 
Solution for the + sign of the dispersion has the form 



=(+) 



= + ^ e^'r L + . /T^ 2 tanh EEM 

v-y^W \ V 2 



e(") = - 



v+^v 2 ~ 4/3 2 
and the density is 



-2/3 2 r 



+ + To) j 

(73) 

v 2 - 4/3 2 tanh + ur + a; )j 

(74) 



1 H 2 seen — -(a; +ur + £0) 



4/3 2 



(75) 
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It shows that velocity of soliton is bounded from below by modulus \v\ > 2(} , 
so the soliton has the "tachionic" character. These results show first that in 
contrast with defocusing NLS with soliton's velocity bounded from above (sub- 
sonic type), the RNLS soliton has velocity bounded from below (supersonic 
type). Another difference is that soliton of defocusing (repulsive) NLS is the 
hole (bubble) like excitation with \ip\ 2 = p < 1, while for the RNLS soliton we 
have the wall like form p > 1. 



7.2 Two soliton solution 

To construct two soliton solution following Hirota [5D] we consider 

g& = a^e" 1 + aje" 2 , h = e" 1 + e" 2 (76) 

Substituting to bilinear equations 

(TD T ± 2kD x + D 2 x ){afe^ + a±e" 2 ) • 1 + 1 ■ (e" 1 + e" 2 ) = 0, (77) 
2{D 2 X + 2/3 2 )(l • (e" 1 + e" 2 )) = 2/3 2 (a+e" 1 + a+e" 2 + a^e^ + a^e" 2 ) (78) 

we have the system 

of {Td T ± 2kd x + d 2 x )e^ + (±d T T 2kd x + d 2 x )e r ' 1 + 
a±(Td T ± 2kd x + d 2 )e" 2 + (±d T T 2kd x + d 2 )e" 2 = 0, (79) 
(kf + 2/? 2 )(e" 1 + e" 2 ) = p 2 [(a+ + + (a+ + a^-] (80) 

Using dispersion 

uf = h ^2k ± y/k?+40^j , (i = 1, 2) (81) 

we have 

+ _ (ujj - 2kh) + kf = +<j) . _ _ (cjy - 2kh) - kf _ J_ _ 
°* (a;* - 2fcfe,) - kf ~ ' ai (fji - 2kh) + kf a+ [ ' 

The last relations imply 

fa - 2kk t ) = kf coth ^ (83) 



and 



so that 



fci =2#)smh— ^ (84) 



{wi - 2kh) = 20 1 sinh ^ (85) 

We note that two signs in dispersion (|8ip correspond in the above formulas to 
the simple replacement 4>i — > First we restrict consideration with the same 
sign for both frequences. In the next order we have the system 

( T D T ± 2kD x + D 2 x )(g± ■ 1+gf • A + 1 • h) = 0, (86) 

(D 2 X + 2/3 2 )(2 ■ h + fx ■ h) = Wtigt + 92 + 9t9l) (87) 
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or for the first equation 

(Td T ± 2kd x + d 2 x )g± + (±d T T 2kd x + d 2 x )f 2 
+{a±[=F(wi - W2) ± 2fc(fci - fc 2 ) + (fci - fc 2 ) 2 ] 
+a^[T(w 2 - wi) ± 2fc(fc 2 - fei) + (ki - k 2 ) 2 }}e m+r ' 2 = 0, (88) 

It implies solution in the form 

gf = a±e" 1+ " 2 , / 2 = 6 12 e'' 1+ " 2 (89) 
Then the second equation of the system implies relations 

a+ = a+a+b 12 = e^+^b^, a^ 2 = a^a^b 12 = e^+^b^ (90) 

so that 

sinh 2 ii^h 

= l^nm < 91 > 

and the first equation is satisfied automatically. As a result we have solution 

e ( ± ) = ±2f , £(l+9t+9t) (g2) 

1 + Jl + J2 

or 

= g*(l + + + b 12 e^+^±(^+^) 

1 + e 1 ? 1 + e'' 2 + 6 :L2 e'? 1 + r ' 2 ^ ' 

For particular parametrization /3o = 1/2, implying ai = 1, we have two soliton 
solution for density p, 

p= 2 , + , ^ ~ (94) 

where 

A ± = sinh 2 ^±^(1 + e m±0 1+e ^±0 2)+ sinh 2 ^LZ^ e ^-H»±Wi+*») j (95) 

Tji = sinh -^jr + [2fcsinh + — sinh^ij]r + rif\ (i — 1, 2) (96) 

If one of the parameters 4n is vanishing or if cf>i — 4> 2 , the solution is reduced to 
the one soliton form. For example if 2 = then 

sinh 2 fy- 

Analyzing two soliton solution in the soliton's moving frames we can see that it 
describes collision of two solitons type (|9"T|) . moving in the same direction with 
the initial position shifts 



1 -"2 
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so that sinh ^ Axi + sinh !f Ax 2 = 0. 

The different form of two soliton solution we obtain if we choose opposite 
signs for frequences (|8Tj) . so that 



u+ = 2/3 (2fcsinhy +/3 o sinh0 1 ), (99) 
uj- = 2/3 (2fc sinh y - p Q sinh X ) (100) 

Then of = e^ 1 , af = e^ 1 , 

a+ = a+a+b 12 = e^~^b 12 , (101) 
«u = &12 = e-* 1+ ^6 12 (102) 

and 

cosh 2 , 

For /3o = 1/2 two soliton solution is 

.D _|_ 

P ~ [ C osh 2 il^h{\ + e m + e r, 2 ) + cosh 2 0l±02 e , )1+ ^]i (J ' J 

where 

B± = cosh 2 (l + e ^+ e *=F*») + cosh 2 ^±^ 2 . e " I +* ± ^^, (105) 

and 

771 = sinh -^-(x — Xx) + [2k sinh ^ + — sinh </>i]t, (106) 
?7 2 = sinh (x - x%) + [2k sinh ^ - i sinh 2 ]t (107) 



It describes collision of two solitons in the form (|97p . moving in opposite direc- 
tion with initial position shifts 



2 cosh - 



1+02 



A*, = -g-ln 1 (i = l,2) (108) 

sinh ^ cosh yi 4 y2 

In Fig.l we show 3D plot of this solution. 
For the velocity field we have 

e (_) e (+) 

u=^-^ = -2* + (109) 
cosh 2 a_ + cosh 2 b_ 

_| 4 4 

cosh 2 + e »'i-0i + e r/2+02) + cosh 2 ^i±^e I 'i+^-0i+02 

cosh 2 a+ + cosh 2 6+ 

Cosh 2 + e r)i+0i + e r )2 -02) + cosh 2 il±^ e 1i+»72+0i-02 



13 



JV 7 



Figure 1: Two Soliton Scattering 

where 

a T = (sinh — e "i=F0i + sinh tl e m±4>^ ( 110 ) 

b T = (Sinh y + Sinh y ) e ^+^T01±02 ( 1U ) 

It has the same phase shift and describe collision of two solitons in the form 
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